Upper critical fields of quasi-low-dimensional superconductors 
with coexisting singlet and triplet pairing interactions 
in parallel magnetic fields 
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Quasi-low-dimensional type II superconductors in parallel magnetic fields are studied when sin- 
glet pairing interactions and relatively weak triplet pairing interactions coexist. Singlet and triplet 
components of order parameter are mixed at high fields, and at the same time an inhomogeneous su- 
perconducting state called a Fulde-Ferrell-Larkin-Ovchinnikov state occurs. As a result, the triplet 
pairing interactions enhance the upper critical field of superconductivity remarkably even at temper- 
atures far above the transition temperature of parallel spin pairing. It is found that the enhancement 
is very large even when the triplet pairing interactions are so weak that a high field phase of par- 
allel spin pairing may not be observed in practice. A possible relvance of the result in organic 
superconductors and a hybrid-ruthenate-cuprate superconductor is discussed. 

PACS numbers: 74.70.Kn, 74.20.Mn, 74.20.-z 



Anisotropic superconductivities have been studied ex- 
tensively in organic, oxide, and heavy fermion supercon- 
ductors. For example, a triplet pairing is confirmed by 
Knight shift measurements in heavy fermion supercon- 
ductors UPt 3 , while NMR experiments suggest a singlet 
pairing with a line node gap in UPd2Al3. On the other 
hand, in a quasi-one-dimensional organic superconductor 
(TMTSF) 2 C10 4 , a line node gap is supported by NMR 
data by Takigawa et al. , while thermal conductivity 
measurements by Belin et al. suggest a full gap super- 
conductivity |||§. In (TMTSF) 2 PF 6 , triplet pairing su- 
perconductivity is supported by recent Knight shift mea- 
surements ||. 

Pairing symmetries can be different in (TMTSF)2C104 
and (TMTSF) 2 PF 6 , but from the similarity of crys- 
tal structures and electronic states in these compounds, 
probably the pairing interactions have the same origin. 
From the phase diagram in pressure and temperature 
plane, d-wave like singlet superconductivity due to pair- 
ing interactions induced by the antiferromagnetic fluctu- 
ations have been discussed HQ. However, it has been 
discussed recently that such pairing interactions contain 
both singlet and triplet channels as attractive interac- 
tions, and even at zero field inter-site Coulomb interac- 
tions might favor a triplet pairing Pairing inter- 
actions induced by antiferromagnetic fluctuations have 
been discussed also in high-T c oxide superconductors for 
the proximity to the antiferromagnetic phase || . 

In this paper, we examine quasi-low-dimensional su- 
perconductors in which singlet and triplet pairing inter- 
actions coexist. In particular, we concentrate on systems 
in which the singlet pairing interactions dominate at zero 
field. We calculate critical fields of superconductivity in 
directions parallel to the highly conductive layers. Be- 
cause of the parallel direction, we assume that our system 
is strongly Pauli limited and the orbital pair-breaking ef- 



fect can be ignored as a first approximation. 

Matsuo and the present author et al. studied this 
problem in a three dimensional system with a spherical 
symmetric Fermi surface in which s-wave pairing interac- 
tions and weaker p-wave pairing interactions coexist [ flO| . 
We found a remarkable enhancement of the critical field 
due to a mixing of order parameters of s-wave and p- 
wave symmetries. The order parameter mixing occurs 
due to appearance of non-zero center-of-mass momentum 
of Cooper pairs stabilized by Zeeman energy. Such an 
inhomogeneous superconducting state is called a Fuldc- 
Ferrell-Larkin-Ovchinnikov state (FFLO or LOFF state). 
It should be noted that the enhancement occurs far above 
a transition temperature of the pure triplet pairing su- 
perconductivity which is estimated in the absence of the 
singlet pairing interactions. We discussed this effect in 
connection with the phase diagram of a heavy fermion 
superconductor JlT| . The order parameter mixing in the 
FFLO state have been pointed out also by Psaltakis et 
al. and Schopohl et al. in an s-wave superconductor and 
in a p-wave superfluid 3 He p^Jl^ . 

In this paper, we extend our previous study to the 
two-dimensional systems and to the anisotropic singlet 
pairing. We assume inter-layer interactions implicitely so 
that the BCS-like mean field approximation is justified, 
while they are weak enough to be neglected in resultant 
mean field equations. In this sense, our systems are quasi- 
two-dimensional. We find that even very weak triplet 
pairing interactions enhance the critical fields remarkably 
also in the present systems. 

The pairing interactions are expanded as 



^(k,k') = -^ffa7a(k)7a(k') 
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where 7 Q (k) are defined by 7d x2 _ B 2 s ) — 
7 Px (k) = k x and so on in cylindrically symmetric sys- 
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terns. We take units with fi — 1 and ks = 1 in this 
paper. We consider two cases: (i) g s > g p > and (ii) 
fjd > fjp > 0, where we have defined g p = g Px = g Vy and 
gd = gd 2 2 — 9d~„- I n each case, the other coupling 

x —y y 

constants are assumed to be zero. The gap function is 
expanded as 



A(k) =^A Q7Q (k). 



(2) 



The gap equations in the vicinity of the second order 
transition are written as 



A Q = g a ^2 KapAp 



with 



(3) 



(4) 



where we have defined £ = h(qcos8 + 1), the angle 8 
between k' and q, q = vpq/2h and h = |/i e H| with the 
electron magnetic moment /i e , It is apparent from the 
above equations that the mixing of the order parameter 
components of odd and even parities occurs only when 
both conditions of q ^ and h ^ are satisfied, which is 
expected from a symmetry consideration in momentum 
and spin spaces. 

In the weak coupling limit, eq.(3) are rewritten in the 
form 
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with T$!} = 2e" / to c /TT • e-Vff-^-C) and 
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where 



cosh y ■ (cosh y + sinh -If ) 



dip 



N a (0) 



(7) 



2tt 



p aa {_0,ip). 



Here, <p is the angle between k' and fc x -axis, and p(0, ip) is 
the angle dependent density of states at the Fermi energy. 
Equations (5) give the upper critical fields h = h(T, q) for 
a given q. The final result of the critial fields is obtained 
by maximizing h(T, q) with respect to q. 
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FIG. 1. The upper critical fields of superconductivity in 
the case (i) in the presence (solid line) and the absence (dot- 
ted line) of the p-wave pairing interactions. T$ /T^ = 0.1 
is assumed. The vertical thin solid line is the transition tem- 
perature of the pure p-wave superconductivity of parallel spin 
pairing. The broken line shows the Pauli paramagnetic limit 
when the FFLO state is ignored. Here, we have defined 
A s0 = 2^e- 1 ^ jv ^°). 

For the case of s-p- mixing (i), we can choose the vector 
q in an arbitrary direction from the symmetry of the 
system. Thus, let us assume the direction of q in k x - 
axis. Then, the p-wave component with A Py ~ k y is not 

mixed with the s-wave component, while A Px ~ k x can 
be mixid. Thus, we only need to calculate the smallest 
eigen value A of the 2x2 matrix 
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where we have defined 



G p = log ■ 



T, 
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(o) g p N p (0) g s N s (0)- 



(8) 



(9) 



The transition temperature is given by T c 

For the case of d-p- mixing (ii), we can use the sym- 
metry to fix the orientation of <i-wave order parameter. 
Thus, let us assume A^ ~ fc 2 — ky. In the absence of p- 
wave pairing interactions, it is known that the optimum 
q is in the direction of fc^-axis (or equivalently fcj,-axis) at 
low temperatures, while it is in the direction of the line 
of k y = ±k x at high temperatures [Jl4|,[l5). In the pres- 
ence of the p-wave pairing interactions, the direction of q 
is not known a priori. Thus, we should assume that the 
two p-wave components (k x and k y ) of order parameter 
can be mixed with the <i-wave components. Therefore, 
we calculate the smallest eigen value A of a 3 x 3 ma- 
trix with elements Mdd, Md Px , M Pa . Px + G p , ■ ■ ■ and so on, 
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which gives the transition temperature by T c = T c y e } 
Here, the parameter G p is defined by 



G p = loj 



T. 



,(0) 
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similarly to eq.(9). 



point of the normal phase, the BCS superconductivity 
(q = 0) and the FFLO state (q ^ 0). It is found that 
the temperature T* increases due to the order parameter 
mixing in the presence of the p-wave interactions from 
the value T* « 0.561 x T c ( ° } in the absence of the p-wave 
interactions. For example, T* 
for T c ( p 0) /T c ( d 0) =0.1. 



0.668 x T c y is estimated 
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FIG. 2. The upper critical fields in the case (ii) in the pres- 
ence (solid line) and the absence (dotted line) of the p-wave 
pairing interactions. Tcp /T^°' = 0.1 is assumed. The verti- 
cal thin solid line is the transition temperature of the pure 
p-wave superconductivity of parallel spin pairing. tp q is the 
angle between q and fc^-axis. At each temperature, the higher 
h on the solid line is the final result of the critical field in the 
presence of the p-wave interactions. Here, we have defined 
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Numerical results are drawn in figures [l] - ^|. Fig 
shows the results for the s-p-mixing, while figures |^ - 
show those for the d-p-mixing. In the both cases, re- 
markable enhancements of the critical field are found at 
temperatures far above the transition temperature of the 
pure p-wave superconductivity of the parallel spin pair- 
ing. In particular, Fig.|| shows that the critical field is 

enhanced considerably even when T$/T£> = 0.01. 

For the d-p-mixing case, the optimum q is in the di- 
rection of k x -axis at low temperatures, while it is along 
the line of k y = ±k x at high temperatures, as shown in 
figures p] - [|. In these figures, the results for tp q = and 
tp q = 7T/4 are drawn by the solid lines. For each direction 
of q, the magnitude q = |q| is optimized to obtain the 
critical field. It is confirmed by numerical calculations 
that q with the other directions are not optimum. 

In Fig.||, the optimum value of q = |q| along the sec- 
ond order transition line is drawn. The temperature T* 
at which q vanishes is the temperture of the tri-critical 
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FIG. 3. The upper critical fields for the case of d-p-mixing. 
= 0.025 is 
the same as in Fig.l 



Tcp /Tty = 0.025 is assumed. The definitions of the lines are 
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FIG. 4. The upper critical fields for the case of d-p-mixing. 
cpV^i? = 0-01 is assumed. The definitions of the lines are 
the same as in Fig.Q. 
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FIG. 5. The magnitude of the optimum wave vector q 
at the upper critical fields for T^'/tJ? = 0.1 when tp q = 
(solid line) and ip q = n/4 (broken line). The dotted lines are 
the result for pure d-wave pairing. 

In conclusion, we examined upper critical field of the 
superconductivity when the singlet and triplet pairing 
interactions coexist. We extended our previous study 
in a spherical symmetric system [ [To| to the quasi-two- 
dimensional systems in parallel magnetic fields and to the 
anisotropic singlet (d-wave) pairing. We obtained similar 
results to our previous results, and find that even very 
weak triplet intereactions enhance the critical fields re- 
markably, and the temperature of the tri-critical point of 
the normal, BCS, and the FFLO states is also enhanced. 
The present phase diagrams coincide with one predicted 
in our previous paper except that the orbital pair- 
breaking effect is not taken into account in the present 
paper. 

In a high-T c superconductor RuSi^GdC^Os, for the 
coexistence of the superconductivity and the ferromag- 
netism (or canted ferromagnetism jl6|), the FFLO state 
has been discussed [jl7],|l8|]. Probably, e?-wave pairing in- 
teractions are dominant in this system, but weak triplet 
pairing interactions must coexist because of the proxim- 
ity to the magnetic phase. Therefore, the present mech- 
anism which enhances the critical fields and T* may play 
a role in stabilizing a bulk superconductivity in this sys- 
tem @. 

In the organic superconductors, the FFLO state has 
been discussed by many authors ||l9|-pT|] to explain the 
high upper critical fields which exceed a conventional 
estimation of Pauli paramagnetic limit (Chandrasekar- 
Clogston limit). As we have discussed above, there is 
a possibility that the antiferromagnetic fluctuations con- 
tribute to the pairing interactions in these systems, and 
they should contain both singlet and triplet channels as 
attractive interactions Q. Therefore, the present mech- 
anism may contribute to stabilizing the superconducting 



phase at high fields. As shown in Fig.[|, the enhancement 
can be very large even when the triplet pairing interac- 
tions are so small that the pure triplet superconductivity 
of parallel spin pairing is not observed in practice. 
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